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We study the transmission and reflection of a plane electromagnetic wave through a two dimen¬ 
sional array of rf-SQUIDs. The basic equations describing the amplitudes of the magnetic field 
and current in the split-ring resonators are developed. These yield in the linear approximation the 
reflection and transmission coefficients. The polarization of the reflected wave is independent of the 
frequency of the incident wave and of its polarization; it is defined only by the orientation of the 
split-ring. The reflection and transmission coefficients have a strong resonance that is determined 
by the parameters of the rf-SQUID; its strength depends essentially on the incident angle. 

PACS numbers: Josephson devices, 85.25.Cp, Metamaterials 81.05.Xj, Microwave radiation receivers and 
detectors, 07.57.Kp 


I. INTRODUCTION 

Recently a layer of metamaterial containing specially 
etched designs, referred to as metasurfaces, was used to 
induce a phase gradient in an incident electromagnetic 
wavejT. This leads to a generalized Snell’s law and the 
control of the transverse structure of the wave front [I]. 
This study was done in the optical domain, it then was 
extended to microwaves by Shalaev et al [2]. The phase 
gradient is due to a plasmon resonance between the elec¬ 
tromagnetic wave and the designs etched on the surface. 
These have to be adapted for each frequency domain and 
are fixed by construction. 

It would be useful to have a system whose response 
could be changed over a significant range of frequencies. 
Such a device exists, it is a split-ring Josephson resonator 
(rf-SQUID) [3-6]. The rf-SQUIDs as basic elements of 
quantum metamaterials were discussed in [THIO]. In the 
experimental study m, it was shown that one can tune 
the resonance frequency of these devices. It should be 
pointed that this system has a discrete energy spectrum 
and a large magnetic momentum m- Then the energy 
of interaction with an external field can be of the order of 
the transition energy between neighboring energy states. 

In this article we show that a film of properly oriented 
rf-SQUID controls the polarization of a wave reflecting 
on the meta-surface. This is similar to a Faraday effect; 
the wave is strongly reoriented at the resonance. We de¬ 
termine the parameters of the reflected and transmitted 


waves in the linear approximation. The reflection and 
transmission coefficients depend on the frequency and on 
the stationary state of the system. In that sense, the 
device is active, its parameters can be modified. 


II. THE MODEL 

We consider a plane wave normally incident on a layer 
of rf-SQUIDs as shown in FigJTJ All the rings in the 
layer are oriented in the same direction given by the 
normal vector n. The interaction of the individual rf- 
SQUID with the electromagnetic field is determined by 
the magnetic flux through the split ring. Therefore the 
orientation of the rf-SQUID controls the parameters of 
the transmitted and reflected waves. This orientation is 
characterized by the angle 6 between the magnetic com¬ 
ponent H and the normal n to the split ring. We assume 
the same homogeneous dielectric layers above and below 
the the film. The dielectric permittivity of this medium 
is 5. The model describing the interaction of electro¬ 
magnetic field with the system of rf-SQUIDs is based on 
Maxwell’s equations and the equation for the response of 
the rf-SQUIDs: 


V x E — — (moH + M) t , 

(1) 

V x H = e 0 eE t , 

(2) 

V • D = 0, V • B = 0. 

(3) 
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FIG. 1: Schematic view of the electromagnetic wave in¬ 

cident on the layer containing split ring Josephson junction 
resonators. One ring is shown, indicating the orientation. 


The magnetization M in the equation 0 is localized in 
the array whose thickness is much smaller than the wave¬ 
length A. Therefore the magnetization can be written as 
follows: 


M(£,r) = a (t)5(z) ~ m (t)n r l5(z), (4) 

a 

where m a (t) are individual rf-SQUIDs magnetizations, 
l « A is the film thickness, n r is the density of rf- 
SQUIDS resonators and m (t) is the magnetic moment of 
the ring |6j. From Maxwell’s equations 0i we obtain 
the wave equation 

H tt - —AH = ——M tt . (5) 

e Mo 

Since the model considered is translationally invariant 
with respect to the x axis, this equation can be presented 
as 

H„ - -H zz = ——M„. (6) 

e Mo 

Let us now consider the Josephson split-ring part. The 
magnetic moment is 


m (t) = Sin. 


(7) 


where I is the current in the loop (11). Combining the 
two equations 0 and 0 we get the final expression for 

M 


M = SIn r l5(z)n. (8) 

Plugging the above expression into the wave equation 0 
yields the wave equation 

H« - -H zz = ~ — SI tt n r l6(z)n. (9) 

£ Mo 

As in m the current in the ring is given by 
I = —L _1 (4> + 27r4>oM), 


where <£> is the flux induced by the electromagnetic field, 
where L is the inductance of the loop, where p is the 
supraconducting phase in the junction and where <f>o = 
^ is the reduced flux quantum. The flux <£> across the 
ring of area S is 


<$> = S H • n 


where n is the normal vector to the ring. This gives us 
the current in the ring 

I = —L~ l (SH • n + 27r$o9?). (11) 


The evolution of the variable p is the same as in m 
except that the current on the right hand side is modified. 
We have 


2^*0 W +2w J[^ +IcSiniP = L (12) 


Plugging (11) into (12) and multiplying by L /( 27t ( I>o ) we 
get 


L Lie 

LCvtt + Jpt + ^ 


sm p = 


27r4>o 


H n — p. (13) 


Our model consists in the wave equation ([9| together 
with the split-ring Josephson equation (13). These can 


be normalized as in m- The natural units of time, flux 
and space are 


cut = 1 /\/LC, 4>o, Zq 


c 


where ujt is the Thompson frequency and zq is the inverse 
of the Thompson wave number. The magnetic field H is 
normalized as 


h = 


US 

27t4> 0 ’ 


r = cj T t, C = z/z 0 . 


In terms of these variables the final equations are 
h rr - h<^ = k (h rr • n + p TT ) 5(C) n. 

p TT + ap T + p + P sin p = — h • n, 
where the parameters <a, /3 and k are 


_ 1 L o LI C 

a r\I C’ 13 27r$o ! 


JS 2 


K = 


HqLzq 


(14) 

(15) 

(16) 

(17) 


The equations ( 15|16 ) are our principal model and we 
now proceed to analyze them. 


III. MICROWAVE SPECTROSCOPY DATA 

As in m we assume that the film is submitted to a 
stationary field that will fix the phase p s . We can then 
write 


( 10 ) 


h = h s + (5h, <p = ip a + Stp 
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where the field dh and the variable 8p are small. Their 
evolution is given by 

Sh rr — 8h^ = k (£h rr • n + 8p rr ) 8(() n. (18) 


Sp jTT + aSp :T + Sp + /3 cos p s 8p = —8 h • n. (19) 

To solve the equations HU we assume the usual 
harmonic dependence 


Sh = e iut i, 6<p = e iut 4>. 


Using these expressions we obtain 


f r = i— M(w)£>(w)(f in • n)n, (33) 

ftr = fin + i^M(w)£>(w)(f in • n)n. (34) 


Equation (33) implies that the polarization of the re¬ 
flected wave is determined by n only. On the other hand 
the direction of the transmitted wave depends on the ori¬ 
entations of n and fi n and on the frequency uj. 

From expression ([33]) the reflection coefficient is 


This yields the following equations 

f££ + cj 2 f = c j 2 k (f • n + (p) 8(() n, 


( 20 ) 


|f r | 2 _ (wk/2) 2 |M(cu) | 2 cos 2 0 
|f in | 2 l + ( W «/2)2|M(a;)|2 ’ 

where 0 is defined by the formula 


(35) 


(j> (uj 2 — uj 2 + aiuo ) = —f • n, (21) 

where we have introduced the resonant frequency uj r 

uj 2 = 1 + /3cosp s . ( 22 ) 

We now set up the scattering experiment by assuming 
an incident field and calculating the reflected and trans¬ 
mitted fields. For £ < 0 we have 

f_(C) = fi„e-^+f r e ia;C . ( 23 ) 


fin • n= |f in |cos0. (36) 

The transmission coefficient is T = I ftr | 2 / |fin | 2 • It is 
^ 1 (wk/2) 2 |M(o;)| 2 c 0 s 2 6 ) 2 

T = ITKflSHF cos e + 


+ i 


.ujk [M(cj) — M*(cj)] 


2 1 + (^/2) 2 |M(cj)| 2 


cos 2 6 = 


= 1-31- 


ImM(c< 


1 + (^/2) 2 |M(cj)| 2 


6. 


(37) 


For ( > 0 the field is 

f+(C) = ftre^. 

At the interface ( = 0 f is continuous so that 


fr = ft, 


We have the following jump condition for 
[fdot = Kw2 ( f tr -n + 4>)n. 


This yields 


fin fr — ftr ZUJK/(ffa ' H T 0)n* 


From the equations (25, 27) it follows 


fin = ftr - i— M(w)](f tr • n)n, 
fr = i^M(w)(f tr • n)n, 


where 


M(co>) = 1 + (uj 2 — uj 2 — iauj) 1 . 


From the first equation we obtain 


(ftr • n) = D(u)( fin • n), 


where 


D(u) 




1 -1 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 



QJ 


FIG. 2: The reflection coefficient IR(cj) from (35) for a plane 
wave normally incident to the film,# = 0. 


Let us analyze numerically the reflection and transmis¬ 
sion coefficients. We choose /3 = 30 so that we are in the 
highly hysteretical case as in [1T5]. Fig. [2] shows the de- 
pendance of the reflection coefficient on uj for a normal 
incidence in the lossless case a = 0, when k, = 1 and 
(3 = 30. The external magnetic field h s is assumed to 
be zero. Then we can take p s = 0 which is the global 
minima of the potential of the equation © (see Fig 2. 
in EH). The reflection coefficient shows a strong reso¬ 
nance for frequencies near coy; this resonance is of the 
Fano type (see EH). The value of the resonance fre¬ 
quency is 

cj r = y/l + /3 cos (f) s ~ 5.56. 
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The metasurface is transparent for small u j; when the fre¬ 
quency is large the incident field is totally reflected. The 
shape of the spectrum does not depend on the specific 
value of (3 and uj r . For instance, the experimentalists 
[T 2 I chose a smaller /3 and obtain a reflection coefficient 
that behaves similarly to the one shown in Fig. [2] 

The magnetic component of the reflected wave is al¬ 
ways oriented in the direction of the normal. On the 
contrary, the polarization angle of the transmitted wave 
f tr depends on the frequency uj and on the orientation 
of the split-ring (angle 9). To analyze this dependance 
assume that the incident field is parallel to the y axis. 



Then the reflected and transmitted fields can be written 
as 


fr = 




From the expression (39) for Ti j2 one obtains the angle 
of polarization rotation of the transmitted wave 'ip in the 
(y, z) plane: 


'ip - arctan 


Re X 2 
Re XI 


(40) 


Fig. [4] shows this angle ip as a function of uj and 9. 





9 ) 


From the relations (30),(32),(33),(34) we get 


Ri—Tx- 1, R 2 — X 2 , 


Ti = D(uj) 


KUJ 9 ^ w ' 

i— sin 2 9M(uj) 


m ikuj sin# cos# . 

X 2 =--- D(oj)M(uj). 


(38) 


(39) 


The transmission coefficient T(u;) is shown in Fig. [ 3 ] 
for three different values of the orientation angle # = 
27r/5, 7 t/4 and 7r/6. 



Cjl) 


FIG. 3: The transmission coefficient 7(uj) from ( |37| ) for a 
plane wave and three values of the orientation angle # . 


Note that the array is transparent 7=1 when uj = uj r . 
Past this frequency, the transmission coefficient asymp¬ 
totically tends to a constant. For a large angle 9 = 27t/5, 
the ring has less influence, since most of the field is trans¬ 
mitted. When the angle is reduced, the incoming mag¬ 
netic field interacts strongly with the rf-SQUIDs so the 
reflected field is stronger. 


FIG. 4: Angle ip(c 0) from (40) as function of uj and 6, a = 0 



9 ) 
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FIG. 5: Angle ip(<jS) from (40) as function of uj and 9, a = 0.05 


Fig- [4] demonstrates the Faraday effect which takes 
place in rf-SQUID metasurface under consideration. The 
polarization rotation angle 'ip depends both on the fre¬ 
quency uj and on the inclination angle of the rf-SQUIDs. 
The angle of rotation of the polarization changes sharply 
near the resonant frequency uj r (gigantic Faraday effect). 
This angle depends strongly on 9] it has a pronounced 
resonance character near uj r . In this case, when 9 changes 
from —7 t/2 to 0 angle ip monotonically decrease from 0 
to —7 t/2 . For 9 = 0 the angle 'ip jumps from —7r/2 to 
7r/2 and monotonically decreases to 0 when 9 changes 
from 0 to 7r/2 (see Fig. [4|. Introducing losses which are 
always present in any practical situation smoothes the 
jump transition, Fig. |5| 
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IV. DETAILED STUDY OF THE RESONANCE 


We study here the characteristics of the resonance of 
the coupled system - rf-SQUID and field. From equations 
(19) and ( [34] ) we can write the time dependent phase 
equation asQBj) 


+ a<p(t) + (1 + /3cos <j> 8 )ip(t) = - f(t ), 
f(t) = 7 ; (y(t) + f(t )) + k cos ujt. (41) 


Recalling the resonant frequency oj r from (22) we replace 
cos ujt i y exp(iu;£). 

so that the system © can be rewritten: The solutions 
of the resulting system have the following form: 


where k = k/2 . Since a and k are small, up to 

second order has following form: 


UJ 


2 

max 


00 r 


or 


2 


— an — 



When ljJ is close to n; r , (u = uj r + A, |A| <0) and a and 
n are small, then 


k 

(j r ^/4A 2 + (a + k ) 2 


(43) 


If K = 0, we recover the classical resonance observed 
for a damped driven linear oscillator, both for 0 and a, 

m ■ 


c p(t ) = Aexp(iut), f{t ) = T>exp(iu;£), 

A = -kC~ \ B = k[(uj?-uj 2 )+iauj\C-\ 

C = id 2 — cj 2 (1 - |a) + iuj [a + f (1 + u 2 - oj 2 )\ . 

Representing A in terms of a phase and an amplitude 
A = aexp(i$) we obtain: 


^/[w2- w 2( 1 -(Ka/2))] 2 +u; 2 [a+( /t / 2 )( 1 +w 2 -w2)] 2 ’ 


tan 0 = 


u> |^a+(/t/2) (1+u; 2 — a; 2 )j 
uj 2 (1 — (kol/2)— 




FIG. 6: Plot of the amplitude a(uS) (left panel) and the phase 
0(uS) (right panel) for three different values of the coupling n = 
0.4, 0.2 and 0.1. The other parameters are uo r — 5, a = 0.01. 


V. CONCLUSION 

We analyzed the interaction of a plane electromag¬ 
netic wave with a two dimensional array of rf-SQUIDs. 
The wave vector of the incident field is assumed orthog¬ 
onal to the array of rf-SQUIDs. All rf-SQUIDs have the 
same inclination with respect to the surface of the array 
and the effective thickness of this ” array film” is much 
smaller then the wavelength. Our main result is that de¬ 
spite this small thickness, the array effectively controls 
the wave reflection and transmission. In particular, it 
changes the polarization of the reflected wave and this 
change is determined only by the orientation of the rf- 
SQUIDs. This effect is identical to the Kerr effect in a 
gyrotropic medium. Here the gyrotropy is introduced by 
the rf-SQUIDs. The polarization of the transmitted wave 
also changes and depends both on the carrier frequency 
and on the inclination angle of the rf-SQUIDs. This is 
similar to a Faraday effect. At the resonance frequency 
we obtain a gigantic Faraday effect. We emphasize that 
the thickness of the array is always much smaller than 
the wave-length. 

This array of rf-SQUIDs acts as a met a-surface that 
controls the polarization of an electromagnetic wave. The 
analysis that we carried out is limited by the linear ap¬ 
proximation. Increasing the incident field will cause a 
larger current in the ring and subsequently a nonlinear 
response of the rf-SQUIDs. We then expect nonlinear 
Kerr and Faraday effects combined with bistability. 


The amplitude a and the phase 6 are plotted as a func¬ 
tion of the frequency uj. Notice how the resonance gets 
sharper for a small coupling n. We also have the typical 
jump in the phase as one crosses the resonance. 

The amplitude a{uj) reaches its maximal value for 


00 max ~ ^2 [^ 2 ( 2cJ r “ “ 2 ) - l] + 


1 

_ 3 k 2 


1 T n 2 [4 T 2(jj 2 — a , 2 — Gna-\- 


k 2 [(uj 2 - 1) [(uj 2 - 1) - 4<n 2 ] + a 4 ]] 


1/2 


(42) 
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